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The phase diagram of the spin-1 bond-alternating XXZ chain is studied numerically. This 
model is appropriate to study the VBS picture and the hidden Z2 x Z2 symmetry concerning 
to the Haldane gap problem. The possible phases are the ferromagnetic, the XY, the Haldane, 
the dimer, and the Neel ones. The critical properties of this model can be interpreted as the 
Ashkin- Teller type reflecting the hidden Z2 x Z2 symmetry. The phase transitions are of the 
Berezinskii-Kosterlitz-Thouless type (XY-Haldane and XY-dimer), of the 2D Gaussian type 
(Haldane-dimer), and of the 2D Ising type (Haldane-Neel and dimer-Neel). We determine the 
phase boundary and universality class using the recent developed techniques by the authors. 
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§1. Introduction 

HaldancB made a fascinating conjecture that the in- 
teger spin 5 Heisenberg antiferromagnetic chain has a 
unique disordered ground state with an energy gap, while 
for a half odd integer 5 it is critical without energy gap 
and belongs to the same universality class as the 5 = 1/2 
case. For 5=1, this-conjecture is practicaUy-.established 
by many numfiricalotiP and experimentao&tiP studies. 
Affleck et. alii proposed the valence bond solid (VBS) 
state for the Haldane gap systems. They studied a 5 = 1 
isotropic spin chain with special biquadratic interactions 
and constructed the exact ground state whose spin cor- 
relation decay exponetially. They also showed the exis- 
tance of a finite gap for the excitation. 

For this problem, den Nijs and Rommelsetl proposed 
the string order parameter, 



ied the 5 — 1 XXZ spin chain with bond alternation, 



string 
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(where a = z, x and (•) means an expectation in the 
ground state) which measures the hidden antiferromag- 
netic order and characterizes the 5=1 Haldane phase. 
By the Goldstone's theorem, if some symmetry is broken 
in the Haldane gap system, it is not expected a contin- 
uous one. Later, the string order parameter was related 
to the hidden discrete Z2 x Z2 symmetry by Kennedy 
and TasakiO. Using a nonlocal unitary transformation, 
they transformed the 5 = 1 spin chains to one model 
which has local interactions and the discrete symmetry 
explicitly. The string order parameters are transformed 
to the ferromagnetic ones, and they relate to the hidden 
Z2 x Z2 symmetry. _ 

In our previous paperM* with Okamoto, we have stud- 



N 



A5J5J +1 ] 



(1.2) 

and found that the phase diagram reflects the hidden 
Z2 x Z2 symmetry. The possible phases are the fer- 
romagnetic, the XY, the Haldane, the dimer, and the 
Neel ones (Fig.l). For 5 — 0, many numerical studies 
have been done. The Haldane-Neel transition occurs at 
A c i = 117 + .^"r^ belongs to the 2D Ising univer- 
sality classBy^BEJ. Near A = A c i, the energy gap 
behaves as |A — A c i|, so the critical exponent v is 1. 
At A = precisely, the Haldane-AF transition occurs 
andi± is of the Berezinskii-Kosterlitz-Thouless (BKT) 
typalJ-*. For this transition, the energy gap closes as 
exp (— const /VA), and in the XY phase the excitation is 
gapless. In the Haldane phase (0 < A < A c i) the hidden 
Z2 x Z2 symmetry is fully broken and the string order 
parameters Cg^ rm g [ol — z i x ) have finite values, while in 
the Neel phase (A cl < A) the discrete symmetry is pa 
tially broken and Cg trmg is non zero but Cg tring = Ol 
Introducing 5, we expect the dimer phase, in which two 
5=1 spins at j = 2i — 1-th and 2i-th sites (5 > 0) 
make a singlet. In this phase, the string order parameters 
are zero, so that the hidden discrete symmetry 
is not broken. We have obtained the phase diagram with 
the same topology of the quantum Ashkin-^lEeller model 
which has the Z 2 x Z 2 symmetry explicitlyLj) . 

But our phase diagram is not consistent with the pre- 
vious one by Singh and GelfancO 1 . They treated the 
model(1.2) by the series expansion and showed that there 
exists a multicritical point of the XY, the Haldane, the 
dimer, and the Neel phases. However, based on the 
Z2 x Z2 symmetry , or the phase diagram of the quantum 
Ashkin- Teller model, there should be a 2D Gaussian crit- 
ical line (Haldane-dimer critical line), which contradicts 
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the phase diagram obtained by Singh and Gelfand. 

In our previous paper we mainly studied the BKT 
transitions (XY-Haldane, XY-dimer). In this paper, we 
investigate the universality class and the critical prop- 
erties of not only the XY-Haldane and the XY-dimer 
transitions but also of the Gaussian transition between 
the Haldane and the dimer phases, and the 2D Ising 
transitions of the Neel and the dimer phases. We clarify 
the relation between the hidden Z2 x Z2 symmetry and 
Haldane-dimer transition, using twisted boundary con- 
ditions. We determine the phase boundary and evaluate 
the critical exponents using the exact diagonalization for 
the finite size systems. 

The organization of this paper is as follows. In the next 
section, we review the quantum Ashkin- Teller model 
which has the Z2 x Z2 symmetry explicitly, and as its ef- 
fective theory we consider the sine-Gordon model which 
descri bes t he information for the critical properties of the 
model(1.2). In section 3, we study the critical properties 
of the XY- Haldane and the XY-dimer (BKT type), the 
Haldane-dimer (2D Gaussian type), and the Haldanc- 
Neel and the dimer-Neel (2D Ising type) transitions using 
the conformal field theory with the numerical analysis. 
The last section is the conclusion. 

§2. Universality class, Z2 x Z2 symmetry and 
U(l) symmetry 

In order to understand the possible phase diagram 
with the Z2 x Z2 symmetry, we review the one dimen- 
sional quantum Ashkin- Teller modelES 

Hat = -^[ojOj+i + t-t? +1 + \<t?<t] +1 t![t? +1 ] 
3 

-/^[aJ + rJ+AaJr/], (2.1) 



where oy ,Tj' are Pauli matrices. This Hamiltonian 
is invariant under a z <-» — a z and r z <-> —t z trans- 
formations respectively (Z% x Z2 symmetry!— For this 
model(2.1), there is a duality transformation!!^ 

i<3 



~x t z h- z ^r z 

r j+l/2 — T j T j+l> T j+ 



1/2 =n 



t<3 



Under this transformation, the Hamiltonian is trans- 
formed as 

H AT (\,f3)=0H AT (\,l/0), 

so at j3 — 1 the system has the self-duality. It is 
known that there is the 2D Gaussian critical line ((3=1, 
— 1/\/2 < A < 1) of continuously varying critical expo- 
nents, and at one end (f3 = 1, A = 1), this line breaks 
up into two 2D Ising critical lines. At another end of the 
Gaussian critical line ((3 = 1, A = — l/v2), it meets two 
BKT critical linesElP , so there is a massless region or the 
"critical fan" in the quantum Ashkin- Teller model. The 
Gaussian critical line separates the fully ordered and the 
fully disordered phases about the Z 2 x Z 2 symmetry, and 
the two 2D Ising critical lines are the boundaries of the 



partially broken phase. ._. 

According to Kohmoto, den Nijs, and Kadanofllia) , the 
quantum Ashkin- Teller model(2.1) can be mapped to the 
5=1/2 XXZ chain with bond-alternation |l|) (A = A, 
S = (1 — (3)/(l+{3)). For this model the Gaussian critical 
line corresponds to the 6 = line [—1/^/2 < A < 1), 
and at the A = 1, 6 = point the bifurcation to the 2D 
Ising critical lines occurs. The universality class of the 
A = 1, 5 = point is_of the level- 1 SU(2) Wess-Zumino- 
Witten(WZW) typeEJ*. 

The effective Hamiltonian of the ID XXZ model 
with bond-alternationE} (or the quantum Ashkin- Teller 
modelEl ) is given by 



H = 



2tt 



dx 



yiv 
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vK(-Klif 



dx cos V2(t>- 
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where n is the momentum density conjugate to </>, 
[(f>(x),U(x )] = iS(x — x ), a is the lattice constant and 
v is the sound velocity. The dual field 8(x) is defined 
as d x 9(x) = 7rn(x). Here we make the identification 
4> = 4>+\/2tt, 8 = 6+^/2tt, which means the U(l) symme- 
try for the 8 field. For the free field theory (y\ = y 2 = 0), 
the scaling dimensions of the operators exp(±im\/2</>) 
and exp(±mv2#) are Km? /2 and n 2 /2K respectively 
(where integer variables m and n are magnetic and elec- 
tric charges in the Coulomb gas pictureES). Since the 
second term of eq.(2.3) is the mass term for the Haldane 
gap systems, we have y\ = for the half odd vakegei-S 
XXZ model, and y x ^ for the integer S XXZ modeO. 

The correspondence between the double sine-Gordon 
model^S) and the quantum Ashkin- Teller model is as 
foUowfl 

1-/3 

yi oc 



1 + /3 



K 



arccos(A) 



Note that the sine-Gordon model(2.3) is invariant under 
the transformation 



7T 

71' 



-> 9, Vi -> -yi, and y 2 -» y 2 , 

(2.4) 

which corresponds to the dual transformation(2.2) of the 
quantum Ashkin- Teller model. 

After the scaling transformation a — > e dl a, we have the 
following renormalization group equations. 



d 1 _ 1 2 1 2 

IlK~ 8 yi + 2 y2 ' 

d yi to K \ 1 

— = (2 - -)y x - -y m , 
^ = (2-2K )m -l yl 



(2.5) 



(This fopn of equations is previously obtained by 
KadanoffEJ.) In Appendix, we show the derivation of 
these equations. Up to the first order of y's, we see that 
yi is an irrelevant field for K > 4 and relevant for K < 4, 
while y 2 is irrelevant for K > 1. Thus in the region 
K > 1, we can neglect the third term in eq.(2.3). In this 
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case, equations (2.5) become 



1 2 
= o2/i> 



dl K 

d Vl fr, K ^ 



which are the recursion relations of KosterlitalEP . There 
isaseparatrix32i^ 1 -81nif- 1 -?^ = 8 + 81n4 (K > 4) 
which separates the infrared unstable region from the in- 
frared stable region, and on this separatrix, BKT tran- 
sition occurs. On this transition point, logarithmic de- 
pendence of finite size spectrum appears, and this makes 
it difficult to extrapolate to the infinite limit in the nu- 
merical calculation. 

The line y± — (1 < K < 4) is the 2D Gaussian 
critical line, which separates two region. For 1 < K < 
4 and y\ ^ 0, y\ flows to infinity. For yi > 0, (4>) 
renormalized to 7t/a/2 as j/i — > +00 and for y-y < 0, 
((f)) —>■ as yi — * —00. 

The equations (2.5) are invariant under the transfor- 
mation yi — > —yi, but not invariant under yi — » — yi. 
For the 5=1/2 bond-alternating XXZ chain(L2), yi > 
and at the SU(2) symmetric point (A = 1, 8 = 0) 
the Gaussian line separates to two 2D Ising critical lines. 
There is another possibility other than the Ashkin- Teller 
multicritical structure, depending on the signature of yi. 
When yi is negative, the Gaussian line connects to a first 
order line (e.g. S — 1/2 XXZ chain with staggered field). 

In order to identify the excitations with those of 
the double sine-Gordon model(2.3), we use the follow- 
ing symmetries. The Hamiltonian(1.2) with the peri- 
odic boundary condition is invariant under spin rotation 
around the z-axis, translation by two-sites (Si — » Sj+2), 
space inversion (Si — > Sn-i+i) and spin reversal (Sf — > 
—Sf, Sj — > —Sf)' Corresponding eigenvalues are 
St = E£=i S f, 9 = *™/ N (n = 0, • • • iV/2 - 1), P = ±1 
and T = ±1 respectively. The corresponding symmetry 
operations in the double sine-Gordon model(2.3) are as 
follows. The operation to the space inversion (P) is 



9 + tt/V2, 



-x, 



(2.6) 



and the operation to the spin reversal (T) is 

</>->-</>, -> -9 + tt/V2. (2.7) 

The correspondence of these symmetry operations is 
summarized in Table 1. 

Assuming conformal invariance, the scaling dimension 
x n is related to the energy gap ofJJae finite size system 
with periodic boundary conditional 



2nv 



(E n (L) - EJL)) , 



(2.8) 



where L is the system size, E g (L) is the ground state 
energy, and v is the velocity of the system. And the 
le^^ig finite size correction of the ground state energy 



E g (L) = e a L--^, 



(2.9) 



where c is the conformal anomaly number. 

In the following section, we calculate these values from 



small size data of the exact diagonalization with the 
Lanczos method. 

§3. Phase transitions and Numerical results 

In this section, wc determine the critical points and 
the universality class. From the hidden Zi x Zi sym- 
metry and the phase structure of the quantum Ashkin- 
Teller model, we expect that the XY-Haldane and 
the XY-dimer transitions are of the BKT type, the 
Haldane-dimer transition is of the Gaussian type, and 
the Haldane-Neel and the dimer-Neel transitions belong 
to the 2D Ising universality class. With the operator 
structure of the sine-Gordon model(2.3) and the renor- 
malization group analysis, we determined the transition 
points numerically. To check the universality class, we 
calculate numerically the conformal anomaly number c 
and several relations between scaling dimensions. 

3.1 Berezinskii-Kosterlitz-Thouless (XY-Haldane and 
XY-dimer) transitions 

First we consider the XY-dimer and the XY-Haldanc 
phase transitions. 

For 6 = case ( X Y-Haldane) , several numerical stud- 
ies have heen done. For the transition point, Botet 

Sd JulianElP estimated A c ~ 0.1, Sakai and Takahashi 
A c = -0.01 ± 0.03, and Yajima and TakahashS 
A c = 0.069 ± 0.003. They expeplexL that the transition 
point is exactly A = (see alsooEO'EEP). In our previ- 
ous paper, we found that the transition point between 
the Haldane and th^YY phases is exactly A = in the 
numerical accuracyElJ . 

The effective Hamiltonian for this transition is de- 
scribed by the sine-Gordon model, 



H 



1 

2^ 



dx 



vK(ttII) 2 + 



dx cos 



K \dx 



(3.1) 



2?ra 2 

Here we neglect the third term of eq.(2.3) which is irrel- 
evant for the BKT transition. For this model, there is a 
special operator, that is the marginal operator 

2 



M = -a 2 K(^Uf + a 2 ^ ( ^ 
K \ ox 



(3.2) 



(this is the Lagrangian density itself for the free field the- 
ory) which changes the scaling dimension of scaling oper- 
ators continuously. For the free field theory, the scaling 
dimension of a/2 cos V2(j) in the second term of eq.(3.1) 
is K/2, and the BKT transition occurs when this dimen- 
sion becomes 2, i.e., the renormalized K is 4. 

Besides the operators M. and \[2 cos V2(f>, v^sin V2(j) 
and exp(=Fi4-\/20) are also marginal at the BKT transi- 
tion point. Recently NomuraEj) studied the structure of 
these operators near the transition point. He found that 
the \[2 cos V2cf> and the marginal operator M. are hy- 
bridized in the course of the renormalization, to satisfy 
the orthogonality condition of scaling operators. We set 
the scaling dimensions of y2 cos \/2(f>- like, y/2smy/2<f>, 
exp(=Fi4\/2#) and marginal-like operators as x\,xi,X3 
and Xq. We denote K and y\ as K = 4 + 2y$ and 
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yi = ±yo (1 + 1) near the BKT critical point (yi = ±j/o)- 
Since the scaling dimension is related to the excitation 
energy as eq.( [2.8] ), Nomura found the finite size behav- 
iors of scaling dimensions up to the first order of yo and 
Vi as 



xi(l) 

x 2 {l) 

x 3 (l) 

x Q (l) 
where we define I 



LAE^L) 

2ttv 
LAE 2 {L) 

2ttv 
LAE 3 (L) 

2ttv 
LAEq(L) 
2ttv 
logL. 



2l/o(0(l 
2/0 (0) 
= 2-y (0. 



= 2 



(3.3) 



yo(0(i+ o*), 



From the recursion relations 
of Kosterlitz(2.6) for t — 0, we have yo(l) = l/log(L/Zn) 
(Lq is a constant), so logarithmic corrections appear. At 
the transition point the scaling dimensions xq and x 3 
cross linearly, and this means that at the critical point 
there exists a degeneracy of excited states. So we can de- 
termine the BKT transition point by the level crossing of 
these excitations, and we can also eliminate the logarith- 
mic corrections by the appropriate average of the excita- 
tions. From table. 1, the first = 0, q = 0, P = T = 1 
[x ] and the Sj, = ±4, q = 0, P = 1 [13] excita- 
tions degenerate on the BKT line. Note that in the 
iSf. = ±1,P = —1 subspace the lowest excitations cor- 
respond to the operators exp(=Fi-\/2#) whose scaling di- 
mension is X4 — 1/8 (77 = 1/4) at the BKT transition 
point. 

In numerical calculations, we study finite rings of N 
sites (N = 8, 10, 12, 14, 16) with periodic boundary con- 
ditions. We determine the velocity v by the current field 
[x = l,q = 2w/L), 

AE(L,k = 2tt/L) 



lim 



2-k/L 



where L is the system size and L = N/2 in the consider- 
ing case. To obtain this value, we extrapolated the value 
of the finite size system as v(L) = i>+const/L 2 . The lead- 
ing finite size corrections of eq. J2.8|) comefcftrn the opera- 



tors L_ 2 i-2l and ((£_ 2 ) 2 + (L_ 2 ) 2 )lBB|3f which are 
the descendant fields of the identity operator. The scal- 
ing dimensions of these operators are 4, and the leading 
correction term of the scaling dimension behaves as 
1/L 2 . (Here we neglect the corrections from third term 
of eq.(2.3) because the scaling dimension of V^cos-v/S 
is 2K, bigger than 4.) In Fig. 2(a) we show the excita- 
tion energies for N = 16, A = —0.5. We determine the 
XY-dimer phase boundary by the crossing point of exci- 
tations SJ = 0, P = T = 1 and S%, = ±4, P = 1. Figure 
2(b) shows the excitation energies for N = 16, A = Q. 
In this case one of the states (Sf. = 0, q = 0, P = T = 1) 
and the eigenstate (Sy = ±4, q = 0, P = 1) are exactly 
degenerate on the whole line, thus in our previous paper 
we have concluded that the XY-Haldane phase bound- 
ary is on exactly A = 0. We determine the multi-critical 
point of the XY, the dimer, and the Haldane phases by 
the crossing point of excitations Sj, — ±4, P — l^a] 
and SJ = 0, P = T = -l[x 2 ] on the A = line. This 



point corresponds to K = 4 and y± — in the continuum 
Hamiltonian (3.1). Figure 3 shows the size dependence 
of this multi-critical point, and we have estimated the 
multi-critical point as A = 0, S = 0.230 ± 0.001. 

In our previous paper, we have presented the size de- 
pendence of (xo + x 2 )/2 and (2x + xi)/3 to verify that 
the transition is of the BKT type. From eqs.(3.3), these 
values eliminate the logarithmic corrections and we have 
seen that these converge to 2 as the system size increases. 
In Fig. 4 we show the extrapolated values along the BKT 
critical line. From eqs.(3.3), the average (3a;o+xi+X2)/5 
also eliminates the logarithmic corrections, where the co- 
efficient 3, 1, 1 is the degeneracy at the critical point. Fig- 
ure 5 shows the size dependence of it, where 1/N 2 be- 
havior comes from corrections of the irrelevant operators 
L_ 2 L- 2 1 and ((L_ 2 ) 2 + (L_ 2 ) 2 )l. In Fig. 6 we show the 
extrapolated values along the BKT critical line. We hav e 
estimated the conformal anomaly number using eq.(2.9), 
as c = 0.998 for A = -0.5, S = 0.583 and c = 1.000 for 
A = 0, S = 0, also considering the finite size correction 
of the x = 4 irrelevant field. Hence we conclude that the 
transitions are of the BKT-type. 

Next, we show the asymptotic behavior near the BKT 
critical point. From eqs.(3.3), we can remove the t linear 
terms as 

l[x (l)+ Xl (l)+x 3 (l)} = 2 + O(y 2 ). (3.4) 



In Fig. 7, we show this value near the critical points 
A = -0.5, 6 = 0.583 and A = 0, S = 0. The linear 
components of t are almost absent. 

Lastly, we consider the XY region where the en- 
ergy gap is zero and the spin correlations decay in the 
power law. The scaling dimensions of the operators 
exp(=FiV2^) [#4] and exp(=Fi2v / 2^) [15] are 24 = l/21f 
and X5 = 2/K, so the ratio of these values is £5/2:4 = 4. 
These operators correspond to the lowest states (Sf. = 
±1, q = 0, P = -1) and (Sf, = ±2, q = 0, P = 1) 
respectively. On the BKT transition point, we have 

^4(0 = I 



1 



yo 

16' 

yo 



Z5(0 — - 4 



(3.5) 



so this ratio is correct up to the first order of the loga- 
rithmic correction y (0 = ^ / \og(L / L n ) . Figure 8 shows 
the extrapolated values of this ratio for A = —0.5. The 
evaluated values are almost 4, except near the BKT 
transition point. The discrepancy near the transition 
point comes from the higher order logarithmic correc- 
tions [0(l/(lnL) 2 )]. 

In the XY region, the product x 2 ■ X3 is 4 up to 
0(l/(lnL) 2 ) (from eqs.(3.3)) near the transition point, 
since the field x 2 is free from the effect of the hybridiza- 
tion. (Contrary to this, the product x\ ■ 23 is affected 
from the hybridization of xq and X\, so there remains 
calO{l/h\L) correction.) In Fig. 9, this product value is 
presented and we can see the expected value. 

In Fig. 10, we show the scaling dimension xq of the 
marginal operator M. of the N = 16 system. On the 
transition point the deviation from 2 is more than 10 %. 
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The deviation of xq from 2 comes from the hybridiza- 
tion with the operator \/2cos^/2(f>. Near the transition 
point we can eliminate this effect of the hybridization 
taking the combination xq + x± — X2 — 2, and so this 
value eliminate the logarithmic correction up to the first 
order of yo. In Fig. 10 we show the extrapolated value of 
Xo(l) + Xi(l) — X2(l), and the discrepancy decreases in a 
few percent. 

3.2 Gaussian (Haldane-dimer) transition 

Next, let us consider the Haldane-dimer transition. In 
the dimer phase, the system has no symmetry breaking, 
whereas in the Haldane phase, the system exhibits the 
hidden Z2 x Zi symmetry breaking. For this transition, 
the operator a/2 cos V2<t> m eq.(2.3) is relevant(4 > K > 
1) and the Gaussian transition occurs at j/i = 0. 
3.2.1 Determination of the transition point 

In the content of the VBS picture, for even site sys- 
tems with periodic boundary conditions, the. Haldane 
gap state of S — 1 systems can be written asES 



N-l 



J] (at&t +1 - 6toJ +1 )|0>, (3.6) 
4=1 

and the dimer state can be 



N/2 



n(4 fe -i4-4-i<4) 2 |o> (<5>o), 



(3.7) 



k=l 



or 



N/2-1 



(^ N bl-bWx) 2 II (4 fc 4+i-4A+i) 2 |o) (*<0), 

fe=l 

where we describe the spin state by the Schwinger 
bosons, that is, at (&t) creates the S — 1/2 | (J.) spin at 
the j-th site. These two states have the same symmetries 
which conserve the Hamiltonian, especially the space in- 
version j — > N — 1 + j (P = 1) and the spin reversal 
a) <-> (T = 1). For the periodic boundary conditions, 
the ground state is singlet in both the Haldane and the 
dimer phases. So we cannot expect the simple finite size 
scaling behaviors near the critical points. 

But the situation is different for another boundary 
condition, i.e., the twisted boundary condition Sff +1 = 
-Sf, S y N+1 = -Sf, S Z N+1 = S{. For this boundary 
condition, the Haldane gap state is written as 

N-l 

(J N b\ + b N a\) (aj&j +1 - &jaj. +1 )|0>, (3.8) 



while the dimer state can be 

N/2 



k=l 



ft 1 

2k-l u 2k 



(3.9) 



or 



N/2-1 



^ N b\+tf N a\) 2 (4 fe 6 



t *t 

2k-\ 



-4«L+i) 2 lo> (s<o). 



k=l 



The Haldane state has P = —1, T = —1, but the dimer 
state has P = 1, T = 1, so these two states have the 
different P and T eigenvalues, and the energy eigenvalues 



cross at the critical point. 

For the sine-Gordon model, we can explain this as fol- 
lows. First, let us consider the following 1-D Hamiltonian 
with periodic boundary condition 

H = H Q + ^- [ dxOi, (3.10) 
2tt Jo 

where Hq is the fixed point Hamiltonian, L is the sys- 
tem length (setting the ultraviolet cutoff as 1), and 
0i is a scaling operator whose scaling dimension is x\ 
(0\ = 0i ). According to Cardyo- 1 , the following size 
dependence of excitation energy up to the first order per- 
turbation is obtained 



AE„ 




(3.11) 



where x n is the scaling dimension of the operator 0„, 
and C n i n is the operator product expansion (OPE) co- 
efficient of the operators 0„ and <D\. We assume that 
the operator 0i is a relevant one (x\ < 2). Up to the 
first order perturbation theory, we find that at the crit- 
ical point A = the scaled gap LAE n does not depend 
on the system size besides finite size corrections of irrel- 
evant operators. But if the OPE coefficient C„i„ is zero, 
the first order perturbation theory is not sufficient and 



we must consider the second order term of A in eq.(B.ll). 
In such a case the scaled gap LAE n has some extremum 
at the critical point. 

In the sine-Gordon model, the relevant operator is 
0i = y2 cos y/2(j), and there is no operator 0„ which has 
nonzero value of C n \ n . It is related with the charge neu- 
trality conditions in the Coulomb gas picture©. (Note 
that the operators e ±% ^l^ is not allowed by the iden- 
tification (j) = (j) + V2tt.) Thus the OPE coefficient in 
cq.( |3.11 ) is zero, and we cannot use the usual finite size 
scaling method. In other words, the Haldane gap sys- 
tems with the periodic boundary condition does not have 
the local "order parameter" which has non-zero value of 
OPE coefficient with the relevant operator in the Hamil- 
toniantf . 

Second, let us see what happens for the twisted bound- 
ary condition Sff +1 ± iS y N+1 = e ±J *(5f ± iSf), S Z N+1 = 
Sf. It has been known that for the primary operator 
exp(iy/2m(j) + iy2n9), the effect of the boundary con- 
dition <f> = 7r is to shift the magnetic charge m by one 
hallH m — ► m — 1/2, n — > n. If the half magnetic 
charges exist, we have the following size dependence of 
the scaling dimensions of the half magnetic charge oper- 
ators v^cos^/v^ [xq] and y/2 sin <j> j y/2 [xt\ asEP 



x 6 (L) 



x 7 (L) 



L(£ 6 (£,$ = 7r)-£ 3 (L,$ = 0)) 
L(E 7 (L,^ = Tr)-E g (L,^> = 0)) 



8 2 U 



2nv 

K/2-2 



(3.12) 



We see that these two energy eigenvalues cross linearly 



6 



Atsuhiro KlTAZAWA and Kiyohide NOMURA 



at yi = 0. 

Note that the sine-Gordon model(2.3) is invariant un- 
der the transformation ( p . 4[ ) . The half magnetic charge 
operator \/2 cos 4>/V% (V2 sin (/>/v / 2)-jCorresponds to the 
operator P = <j z t z (P = a z f z )c3 of the quantum 
Ashkin- Teller model. Under this transformation, the half 
magnetic charge operators \/2 cos 4>/^/2 and \[2 sin 0/V2 
are transformed as 



V^cos -^=c 

\Fi sin —=c 
y/2 



—V2sm — = 
V2 cos —^<t> 



(3.13) 



Thus at the point y± = 0, the system has the self-duality 
in the quantum Ashkin- Teller language. By the invari- 
ance of eq.( [2.4| ), the third term of eq(2.3) does not affect 
the determination of the crossing point for eq.(3.12). 

Numerically, with $ = n boundary condition we de- 
termine the phase boundary of the Haldane and dimcr 
phases for N — 8,10,12,14,16 systems. In Fig. 11, wc 
show two low lying energies of the subspace & z = 
with the boundary condition S% +1 = —Sf, S%+i = 
—Sf, S^ +1 = Sf for A = 0, 0.5 and 1, which corre- 
spond to E 6 (L,§ = 7r) - E g (L, $ = 0) and E 7 (L, $ = 
7r) — E g (L, <f> = 0). From figures of A = 0.5, 1.0, we see 
that in the Haldane phase P = T = — 1 state is lower 
than P = T = 1 state but in the dimer phase the relation 
reverses as expected. 

However, contrary to the S — 1/2 case, since the self- 
duality at y\ = is not exact for 5=1 bond-alternating 
system, the crossing point of E§ and E 7 has a size depen- 
dence as 0(N~ 2 ) (or 0(N~ 4 )) (see appendix B). Figure 
3(A = 0) and Figure 12(A = 0.5) show the size depen- 
dence of the crossing points. In fact, assuming the size 
dependence of the crossing point S C (N) = 5 c (oo) + aN~ b , 
we can check the power b, using the relation 



In 



S c (N+2)-S c (N) 
S c (N)-S c (N-2) 



hi 



( N+ 1 
\ N-l 



" 1 h 



(3.14) 



Figure. 13 shows the exponet evaluated with this method. 
Hereafter we extrapolate the critical point as S C (N) — 
S c (oo) +aN- 2 + bN- i . 
3.2.2 Check of the universality class 

In Fig. 14, we show the value K estimated from x\,x 2 , 
X4, and xq of the N = 16 system. These values are 
almost consistent except near A = 1 where the effect of 
the logarithmic corrections of the level-1 SU(2) WZW 
modclEU appears. 

We estimated the conformal anomaly number at the 
critical point A = 0.5, 5 = 0.2524 as c = 0.998. In 
Table 2, we show some extrapolated scaling dimensions 
for A = 0.5. The estimation of K is consistent in each 
others. 

For the Gaussian transitions the leading finite size cor- 
rections of eq.(|2.8j) come from the operators L_ 2 £-2l 
and ((i-2) 2 + (L-2) 2 )l or \/2 cos \/8cf), of which the 
third operator is included in the third term of eq.(2.3). 
The scaling dimensions of these operators are 4 and 2K 
(1 < K < 4) respectively. By the third term of eq.(2.3), 



at y± = only the scaling dimensions x\ and xi are 
corrected as 



K 






~2 


+ fi 




K 






~2 


-fi 





2K-2 



2K-2 



We can eliminate this correction by summing these two 
dimensions. To improve the precision of the 1-parameter 
scaling, hereafter we try to eliminate the effect of the yi, 
y 2 fields, at least first order. 

We show the extrapolated values of Xq on the Haldane- 
dimer transition lines in Fig. 15. The discrepancy is 
about 30 % at A = 1. Figure 16 shows the extrapo- 
lated values of 



X 4 ■ (X 6 + X 7 ) 



1 

2K 



K K x 



and the obtained values are consistent with 1/8 except 
near A = 1, where the discrepancy (about 15 %) comes 
from the logarithmic corrections of higher orders. In 
Fig. 17, we present the ratio X4 and x§. Figure 18 shows 
the extrapolated values of (xi + X2) ■ X4 which should be 
4. The summation of X\ and x 2 eliminate the correction 
from the third term of eq.(2.3). In Fig. 19, we present the 
ratio (x\ +X2)/ (xe+x 7 ) where the numerator is the value 
with the periodic boundary condition and the denomi- 
nator is with the $ = tt twisted boundary condition. 

Last of this subsection, we remark that the method to 
determine the transition point can also apply to deter- 
mine the Gaussian fixed point (yi = 0) in the XY phase. 
In Fig.l, wc also show the Gaussian fixed line. We find 
that the line ends at the point A = — 1, 8 = 0. Near this 
point, the Gaussian line behaves as S cx (1 + A) 1 / 4 (see 
Fig.20). 

3.3 Isotropic case 

The isotropic case A = 1 has been studied in rela- 
tion with., the Affleck-Haldane prediction. Affleck and 
HaldaneH Elf considered the A = 1 case for the arbi- 
trary S case, mapping onto the nonlinear 0(3) er-model. 
They showed that the topological angle 9 is given by 
8 = 2irS (1 — S) and the system should be massless when 
8/2ir is half odd integer. Numerically several authors es- 
timated the transition point S c . For S p-i-case, they 
concluded A, ~ 0.25 ± 0.01 and c = lBBEP. Tot- 
suka et. alS3 extensively studied this case, and obtained 
the critical dimensions up to the logarithmic corrections 
and identified the universality class as the SU (2) level-1 
Wess-Zumino-Witten model. 

For yi = 0, the RG equations are given by 

d_ J_ 

JlK 
dy2 
dl 

yo/2 near K 

dyp 
dl 



1 2 

2%, 



2(1 



Denoting K = 1 



K)y 2 . 

1, we have 



(3.15) 
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dy 2 
dl 



The WZW^oint is y = 1 
and ScrmlzefP , and Nomura 
scaling dimensions 



According to Giamarchi 
, we obtain the following 



x (l) = 2-yo(l)[l + -t 



xi(l) 
x 2 (l) 



1 3 ,„ / 2 

- 2 + - 4 yo(D i + 3* 



1 1 



2 4 

2:4(0 = 

15(0 = 2-^,(0 



yo(0(i + 2t), 



(3.16) 



1 1 m 

2 - 4^0(0) 



x 8 (0 = 2 + 2y o (0 ^1+3* 

aj 9 (0 = 2 + ito(0, 

where i is defined as y 2 = yo (1+t), and we define the scal- 
ing dimensions of y2 cos v8</> and v2 sin v80 as xg and 
£9 respectively. At t = we have j/o(0 = l/l°g(L/Lo) 
and the eqs.(3.16) are consistent with the non-Abelian 
bosonization theory^ . 

From these equations, we see that (x\(l) + 3xi(l))/A, 
(x (l) + x g (l))/2, (2x (l) + x 8 (0)/3, and {5x + 3x 9 + 
x$)/9 do not depend on the logarithmic corrections at 
the Wess-Zumino-Witten point. 

We obtained the transition point with the same 
method in the previous subsection. The obtained value 
is S c — 0.2598. The conformal anomaly number is 
estimated as c — 0.990. The extrapolated value of 
(xi(0 + 3x4(0)/4 is 0.495. and the extrapolated values of 
(aj (0+a;9(Z))/2, (2s (0+*8(0)/3. and (5x +3a; 9 +a; 8 )/9 
are 1.982, 2.026, and 1.997 respectively. 

3.4 Ising (Haldane-Neel and dimer-Neel) transition 

Lastly, we study the Neel-Haldane and the Neel-dimer 
transitions. Because in the Neel phase, the hidden Z 2 x 
Z 2 symmetry is partly breaking, we see that the broken 
symmetry of these transitions is the Z 2 symmetry. 

On the Neel-Haldane transition (A = 1. 17 ± 0.02, S = 
0), Nomurali-$ calculated the critical exponents as v — 
0.98 ± 0.007, = 0.126 ± 0.007,?? = 0.253 ± 0.002, with 
a (large cluster decomposition) Monte Carlo method. 
Later Sakai and TakahashiH reexamined this transi- 
tion using the phenomenological renormalization group 
method and finite-size scaling technique, and obtained 
v = 1.02 ± 0.05, 7/V = 1.76 ± 0.01, and n = 0.23 ± 0.01. 
These results are consistent with the 2D Ising univer- 
sality class (v = 1, /3 = 1/8, ry = 1/4, 7 = 7/4). So 
here we consider the Neel-dimer transition and check its 
universality class. 

In c < 1 conformal field theory, the current field does 
not exist, so we cannot determine the velocity by the 
same method as the c = 1 case. 

Assuming the 2D Ising transition, we expect that the 
first excitation state has the Sj, = 0, q = 0, P = T = — 1 
symmetry, and this corresponds to the order parameter a 



of the Ising model, whose scaling dimension is 1/8. The 
scaling dimension of the level- 1 descendant field L_ier is 
1 + 1/8, so we can determine the velocity of the system 
by 

Using this velocity we can determine the conformal 
anomaly number c and scaling dimensions of several ex- 
citations. 

Numerically we determined the phase boundary by 
the phenomenological renormalization group method, 
that is, evaluating the crossing point of the scaled gaps 
NAE(N) and (N + 2)AE(N + 2) for N = 6, 8, 10, 12, 14 
systems with periodic boundary conditions. Using the 
velocity ([U7j), 

we estimate the conformal anomaly num- 
ber as c = 0.4997 and the scaling dimensions as x a — 
0.126 and x t = 1.001, for A = 2.0, 5 = 0.683. In 
figure 21, we show low lying scaling dimensions of the 
dimer-Neel transition point, where circle shows Sij, = 0, 
P = T = 1, and cross shows Sf. = 0, P = T = — 1 states. 
These values are consistent with the prediction from the 
Ising universality. 

§4. Conclusion 

We discussed the critical properties of the S — 1 bond- 
alternating XXZ spin chains. The expected phases are 
the ferromagnetic, the XY, the Haldane gap, the dimer, 
the Neel phases. First, we considered the universality 
class based on the hidden Z 2 x Z 2 symmetry. Then we 
have taken the analogy from the quantum Ashkin- Teller 
model which has a Z 2 x Z 2 symmetry explicitly. The ef- 
fective model of the quantum Ashkin- Teller model is the 
double sine-Gordon model, and using the information of 
it, we determined the XY-Haldane, the XY-dimer and 
the Haldane-dimer phase boundaries numerically. For 
the XY-Haldane and the XY-dimer (BKT) transition, 
the renormalization group aspect is important and it is 
difficult to apply the simple finite size scaling method due 
to the logarithmic corrections. While for the Gaussian 
transition between the Haldane and the dimer phases, 
we adapted the twisted boundary condition to obtain the 
preferable operator structure and determined the transi- 
tion point by a level crossing. This method can also be 
applied to determine the Gaussian fixed line in the XY 
phase. We found that the hidden discrete symmetry is 
crucial for the phase transition and the topology of the 
phase diagram. We also identified the universality class 
numerically using the conformal field theory and elimi- 
nating the correction of the (marginal) irrelevant fields. 

It is interesting to consider the arbitrary spin case. Os- 
hikawacS studied the arbitrary S Heisenberg chain with 
bond-alternation in the VBS picture. It was pointed out 
that the successive dimerization transition occurs in the 
Affleck- Haldane prediction^ . Oshikawa calculated the 
string order parameter, and concluded that for the suc- 
cessive dimerization, the breakdown of the hidden Z 2 x 
Z 2 symmetry occurs. Considering the phase diagram 
of the quantum Ashkin- Teller and the Haldane's conjec- 
ture, we predicted in our previous paper that for arbi- 
trary S XXZ spin chains with bond-alternation, there 
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are 25 + 1 BKT lines, 25 2D Gaussian lines, and 25+1 
2D Ising lines in the region — 1 < 5 < 1 (summarized in 
table 3). For <i = 3/2 isotropic case A = 1, Yajima 
and TakahashiEj) studied with the density matrix renor- 
malization group method evaluate the transition point 
as 6 = 0,±0.42. 
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Appendix A: Renormalization group equations 

We derive the renormalization group equations(2.5). 
Let us consider the following Euclidean action, 

where So is a fixed point action, a is a ultraviolet cutoff, 
r = [vt, x), and z — vt + ix, z = vt — ix,. . We set the 
scaling operator O a as the normalized one as 

(0„(^i)<V*»>sO>o = ^rr 5Z-, 

f zi-z 2 Y a ( zi-z 2 \ z " 

(A-2) 

where Z\ = vt + ix, Z\ = vti — ix\ and h and h are 
the conformal weights and the scaling dimension of O a 
is x a h a -\- h a . 

Accpijding to ZamolodchikovcJ , and Ludwig and 
CardyeS), we have the following one loop renormal- 
ization group equations for the scaling transformation 



a = e dl a 



dX 
dl 



= (2- x a )\ a -J2 6 K 



,K ~ ^ lA /3A 7 , (A-3) 



where 



h a — hp — h~, and C Q ^ 7 is the operator 
product expansion coefficient of O a , Op, and 7 . 
For the double sine-Gordon model(2.3), we define 

2 / n / \ 2" 



So 



1 



2irK 



d 2 r 



d<j> 
vdr 



Ox 



Oo = 



K 



vdr 



dx 



Ox = v^cosv^, 

2 = v^cosV^, 

where Oq is proportional to the Lagrangian density. Us- 
ing 

(4>(z, z)4>(0, 0)) = -y log M, 



and Wick's theorem, we can easily derive the following 
operator product expansions 

z)O (0, 0) = -^^-L-OtiQ, 0) + ■ • ■ , (A-4) 



2 (z, z)O (0, 0) = -2K—^O 2 (0, 0) + ■ • • , (A-5) 
\z/a\ 2 



1 (z, zWO, 0) = - 1 lz/a * K/2 _ 2 Oo(0, 0) (A-6) 
1 1 



V2 \z/a\ 



^O 2 (0,0) + 



where '• • •' dose not include Oo, 0\, and 2 . From these 
equations we have 

K 

Cm = Cioi = Cuo = — 

Cq22 = C202 = C220 = — 2if, 



C112 — C121 — C211 — —7=. 



Using these OPEs and eq. ( |A-3| ) , we can derive the renor- 
malization group equations(2.5). 

Appendix B: Correction of the Gaussian fixed 
line 

Our method in §3.2 to determine the Gaussian fixed 
line is affected by the descendant fields (critical dimen- 
sion x — K/2 + 2n; n:integer) of the \/2cos\/20 field, 
since generally the self-duality (in the Ashkin- Teller lan- 
guage) is not exact for the finite size system. Con- 
sidering correction terms from the descendant fields of 
V2cosV20 and L 2 L 2 1, ((L 2 ) 2 + (Z 2 ) 2 )l (x = 4) irrele- 
vant fields, we obtain 



(B-l) 



/o \ K/2-2 oc /o \ K/2+2n-2 

^ ' n—l 



since the OPE coefficient 

(\/2cos (j)/V2\V2 cos V20|\/2cos (f>/V2) 



changes sign under the transformation ( |2.4| ) . Therefore, 
the crossing points behave 

-2n 



yf 



(B-2) 
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Fig. 1. Phase diagram in the A — <5 plane. The XY-dimer and the 
XY-Haldanc phase boundaries are of the BKT type, the dimer- 
Haldanc boundary is of the 2D Gaussian type, and the Neel 
phase boundaries arc of the 2D Ising type. We also show the 
Gaussian fixed line (j/i = 0) in the XY region (x). 
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Fig. 2. Excitation energies for (a) N = 16, A = —0.5 and (b) 
N = 16, A = 0. O's arc 5f, = 0, g = 0, P = T = 1 excitations, 
O is Sj,, q = 0, P = T = — 1 excitations, and + is the SJ, = ±4, 
q = 0, P = 1 excitation. 
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Fig. 4. The extrapolated values of (2xq + rri)/3(0) and (rro + 
a; 2)/2(0) along the BKT transition lines. We also show the 
bare values of xq(x), aii(x), and X2(0) of the N = 16 system. 
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Fig. 3. The size dependence of the BKT multicritical point. O: 
the crossing point of (Sf, = 0, k = 0, P = T = -1) and (Sf, = 
±4, k = 0, P = 1) excitations; Circles arc obtained by the system 
with the twisted boundary condition(§3.2). 



Fig. 5. The size dependence of (3xo +xi + X2)/5 as a function of 
l/A r2 .This behavior comes from the x = 4 irrelevant operators. 
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Fig. 6. The extrapolated values of (3:ro + x\ + X2)/5 (O) along 
the BKT transition lines. 
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Fig. 9. The extrapolated values of X2 ■ £3 in the XY phase. 



Fig. 7. The extrapolated values of (xo + xi + X3)/3 near the BKT 
critical point(O). We also show the bare values of xo(x), a;i(D), 
and 2:3 (+) of the N = 16 system. 
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Fig. 8. The extrapolated values of £5/2:4 in the XY phase. The 
transition point is A = —0.5, <5 = 0.583. 
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Fig. 10. +: The scaling dimension of the marginal operator of 
the A = —0.5,N = 16 system. O: The extrapolated values of 
%0 + xi — X2 which eliminates the logarithmic corrections. 
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Fig. 13. The evaluated value of the exponent b form eq.(25) for 
N = 14. 



Fig. 11. Energy difference E 6 (it) - E g (0) (O) and E 7 (w) - 
E 9(0)(O)- O is the state with ^ S z = 0, P = T = -1, and 
O is the state with J2 S z = 0, P = T = 1. X 
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Fig. 14. The obtained values of K for N = 16 systems from 
xi(0), 2:2(0)1 Z4(rj), and xa(+). 
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Fig. 12. Size dependence of the crossing points. The extrapo- 
lated value is 8 C = 0.2524. 
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Fig. 15. The extrapolated values of xq which is the scaling di- 
mension of the marginal operator. 
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Fig. 16. The extrapolated values of X4 ■ (x$ + x 7 ) on the Haldanc- 
dimer critical line. 
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Fig. 17. The extrapolated ratio between X5 and 24. 
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Fig. 18. The extrapolated values of (x\ + X2) ■ X4. 
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Fig. 19. The extrapolated values of (x\ + X2)/(x% + x-j). 
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Fig. 20. The Gaussian fixed line in the 5 4 — A plane. K takes a 
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Fig. 21. The low lying scaling dimensions of the dimer-Neel tran- 
sition point A = 2.0, S = 0.683. Horizontal line is the momentum 
xL/2-k. O is the scaling dimensions of P = T = 1 states, and 
X is of P = T = — 1 states. 
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x 3 
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exp =F«2 y/20 
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X7 
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V2 cos 2^2(t> 


x s 
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V2sin2%/20 


Xg 



Table 1: Correspondence of symmetry operations, 
defined with the $ = tt twisted boundary condition. 



:These two operators are 



scaling dimension 


21,0 = 1/2K 


rco.i = K/2 


^0,1/2 = K/8 




0.179 


1.41 


0.350 



Table 2: Scaling dimensions at the critical point A = 0.5, S = 0.252. Here we 
extrapolated the corrections from the irrelevant field L-iL^X (x = 4). 
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2-D Gaussian line 


BKT line 
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Table 3: The number of transition lines for arbitrary spins. 
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